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Abstract 

Using the atmospheric neutrinos to probe the density profile of the Earth depends 
on knowing the angular distribution of the neutrinos at production and the neutrino 
cross section. This paper reviews the essential features of the angular distribution 
with emphasis on the relative contributions of pions, kaons and charm. 

1 Introduction 

In principle atmospheric neutrinos can be used to measure the density profile 
of the Earth (integrated along a chord) by comparing the observed rate 
of upward-moving neutrino-induced muons as a function of zenith angle 
with what is produced in the atmosphere. The measurement depends on 
several factors: the zenith angle dependence of neutrino production in the 
atmosphere, the differential neutrino cross section as a function of energy, 
properties of muon propagation in the medium surrounding the detector and 
the angular and energy resolution of the detector. Gonzalez-Garcia, Halzen 
and Maltoni [2008] show that it may be possible to measure some aspects of 
the core-mantle transition region with the IceCube detector currently under 
construction at the South Pole [Karle et ah, 2007]. A data accumulation of 
order ten years or more would be needed. 

The critical energy range for such a measurement is 10 < E v < 100 TeV 
where the absorption of the Earth becomes important, first for neutrinos 
coming straight up through the diameter of the Earth, and at higher energy 
for larger angles. Hoshina [2008] discusses the angular resolution of IceCube 
for neutrino-induced muons from muon neutrinos in this energy band, where 
more than 50% of events can be reconstructed within one degree of their 
true direction. The charged current neutrino cross sections in this energy 
range are known to within ±2% [Cooper-Sarkar and Sarkar, 2008]. In this 
paper, I focus on the angular distributions of muon neutrinos produced in 
the atmosphere and their uncertainties. 
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2 Angular dependence of atmospheric neutrinos 



When cosmic ray protons and nuclei enter the atmosphere they interact and 
produce cascades of secondaries, including charged hadrons, some of which 
decay to produce muons and neutrinos. Muon neutrinos of high energy are 
produced primarily in the decay chain 

K± -> m ± + 

\ (1) 
e ± + v e (V e ) +F M (f M ), 

There are similar expressions for decay chains initiated by pions and charmed 
hadrons. In the energy range of interest here, decay of muons is rare and 
its contribution to muon neutrinos can be neglected. Moreover, the contri- 
bution from kaons is large (~ 80%) as compared to the contribution from 
decay of charged pions (~ 20%). The contribution from charm decay is 
small, but may become significant in the energy region of interest here, as 
discussed below. 

Neutrino fluxes are determined by a combination of four factors: 

1. the energy spectrum and composition of the primary cosmic radiation, 

2. production of pions, kaons and charmed hadrons in collisions of nucle- 
ons and other hadrons with nuclei of the atmosphere, 

3. kinematics of hadron decays to neutrinos, and 

4. the density profile of the atmosphere. 

The primary spectrum can be approximated by a sequence of power laws. 
The energy range relevant for 10-100 TeV neutrinos is of order 30 TeV to 
3 PeV, where the spectrum of nucleons is approximately 

dN (2) 



with 7 « 1.7 [Gaisser & Stanev, 2008]. 

For each of the main neutrino parents (i), there is a critical energy, 

€i = niiC 2 x (-^) , (3) 
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where ho is the scale height of an exponential approximation to the atmo- 
sphere and mj and r, are respectively the mass and rest lifetime of a particle 
of type i = vr ± , K ± , D ± . When 

E v > sec(6>) x ei (4) 

re-interaction of the parent hadron is favored over its decay. Approximate 
numerical values are given in Table 1. 

The energy dependence of the competition between decay and re-interaction 
of the different parent mesons affects the neutrino energy spectrum in an 
essential way. These effects can be displayed explicitly using analytic approx- 
imations for the neutrino energy spectrum. In the power-law approximation 
and at high energy, the spectrum of neutrinos, (f) u , is approximated well by 



ME„) = 4>n{e v ) 



Any 



1 + B nu cos(6)E u /e 7T 1 + B Kv cos(9)E u /eK 

^charm v \ 
1 + -Bcharm v COs(9)E u / e c harm J 

where 4>n{E u ) = dN/d\n(E u ) is the primary spectrum of nucleons (N) 
evaluated at the energy of the neutrino [Gaisser, 1990]. The factors -A.^ 
contain the physics of meson production weighted by the spectrum and the 
decay kinematics contained in Equation [TJ 
As an example, 

Z Nw (1 - r 7r )T 



1 - Znn 7 + 1 
Here 



Oa J dx 

is the spectrum weighted moment for the interaction process in which a par- 
ticle a interacts with a nucleus in the atmosphere and produces a secondary 
particle b that has a fraction x of the lab energy of the projectile. The fac- 
tor involving r v = m^/m^ is the spectrum weighted kinematic factor for the 
decay ir — >• fi + v and 7 ~ 1.7 is the integral spectral index in a power-law 
approximation to the primary cosmic ray spectrum. 

At low energy, the neutrino spectrum follows the same power law as the 
parent hadrons. At high energy, the probability of decay of a charged pion or 
kaon is suppressed relative to hadronic interaction by a factor proportional to 
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Table 1: Critical energies in GeV 



^7r ^charm 

115 850 ~ 5 x 10 



the energy of the meson. This suppression is represented by the denominator 
of each term in Eq. [5j Asymptotically at high energy the neutrino spectrum 
has an extra factor of 1/E relative to the spectrum at low energy, which 
is proportional to the primary cosmic-ray spectrum. This steepening of 
the neutrino spectrum occurs at higher energy for larger zenith angles as a 
consequence of the sec(9) factor, which occurs because particle production 
occurs higher in less dense atmosphere for large angles. The factors Bi v 
in the denominators are ratios of spectrum-weighted kinematic factors for 
meson decay to neutrinos to account for the fact that the decay occurs from 
a steeper spectrum at higher energy. Close to the horizontal (9 > 70°) the 
curvature of the Earth is significant and the secant is to be evaluated at the 
local zenith angle where the particle production occurs, which is less than 
the zenith angle of the trajectory at the detector [Lipari, 1993]. 

Parameterizations of the form of Eq. [5] represent detailed numerical and 
Monte Carlo calculations rather well. The corresponding parametric equa- 
tions for muons can be compared directly to the many measurements of 
atmospheric muons. Such a comparison is made in Fig. 24.4 Gaisser & 
Stanev [2008], for example. (See also Fig. 4 of Gaisser [2005]). 

Uncertainties in the angular dependence arise primarily from uncertain- 
ties in the level of kaon production for E v ~ 10 TeV and at higher energy 
from the more uncertain level of charm production. Ideally one would simply 
measure the downward atmospheric neutrino flux and compare it with the 
upward neutrino flux from the opposite direction and use the ratio Earth-in 
to Earth-out to measure the density profile in a way that is independent of 
the intrinsic angular dependence of the atmospheric neutrino beam. In real- 
ity this unfortunately cannot be done because the relatively high intensity of 
downward atmospheric muons hides the downward atmospheric neutrinos. 
Because of the close genetic relation between and n of Eq. [1] and the 
corresponding equation for pions, a measurement of downward atmospheric 
muons provides some control of the angular dependence. This constraint, 
while valid and worth pursuing, is of limited practical use in this context be- 
cause of the fact that most muons come from decay of charged pions while, 
in the high energy range of interest here, most neutrinos come from decay 
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H. + p. and + flux from pions and kaons 
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Figure 1: Fraction of muons and muon neutrinos from pion decay and from 
kaon decay vs. neutrino energy. Solid lines for vertical, dashed lines for 
zenith angle 60°. 

of kaons. 

3 Neutrinos from kaons 

Each term in Eq. [5] exhibits the same characteristic dependence on energy 
and zenith angle. Kaons become increasingly important with increasing 
energy because ex > e?r- This means that the contribution from kaons 
retains a spectrum close to that of the primary cosmic rays while the pion 
contribution is already beginning to steepen. The location of this transition 
depends on the K/ir ratio at production and is of the order of 100 GeV, as 
shown in Fig. 1. It occurs at higher energy for larger zenith angle. Note 
that kaons never become the dominant source of atmospheric muons. This 
is largely a consequence of the fact that the muon mass is not much less 
than its parent pion so the carries only a small fraction of the energy in 

Uncertainties in the calculated intensity of atmospheric neutrinos at high 
energy due to uncertain knowledge of hadronic interactions were estimated 
by Agrawal et al. [1996]. The single largest source of uncertainty in the TeV 
range and above is from the uncertainty in kaon production, in particular in 
the factor Z pK +. This uncertainty was estimated as ±13% in the 10 - 100 




5 



TeV range. 

The recent measurement of the /i + / yT ratio above 1 TeV with the 
MINOS far detector [Adamson et al, 2007] reflects the importance of kaon 
production in an interesting way. The charge ratio increases as kaons become 
dominant, from which one can infer that the K + / K~ ratio is larger than the 
7r + / 7r~ ratio. It is likely that the specific process of associated production 

p + air — > K + + A + anything (6) 

plays an important role here. Adamson et al. have shown how their mea- 
surement can constrain the parameters that describe the production of kaons 
relative to pions and their charge ratios. Because of the genetic relation be- 
tween atmospheric fi and z/u, the MINOS measurement, in combination with 
measurements of K / tt production at high energy, has the potential also to 
make the prediction of the angular distribution of atmospheric muon neu- 
trinos more precise. 

4 Neutrinos from charm 

Although the charm channel for atmospheric neutrinos is formally the same 
as for pions and kaons, there are two large quantitative differences. One is 
that production of charm in hadronic interactions is very much smaller than 
production of pions and kaons. The other is that the critical energy is so 
high (because of the short charm lifetime) that muons and neutrinos from 
decay of charmed hadrons will continue in the low-energy regime of Eq. [5] 
with the same, relatively hard, spectrum as the primary cosmic-ray nucleons 
and with an isotropic angular distribution. 

Because of its harder spectrum, the charmed contribution will eventu- 
ally become the dominant source of atmospheric v^, adding an isotropic 
"prompt" component to the secant (6) dependence of the kaon and pion 
contributions. The transition energy is at some energy E v > lOTeV, which 
therefore may be in the energy range relevant for Earth tomography with 
neutrinos. Because the level of charm production is so uncertain, this is a sig- 
nificant source of uncertainty in the angular distribution of the atmospheric 
neutrino beam in this high energy range. Gelmini, Gondolo & Varieschi 
[2003] assemble constraints from measurements of the angular and energy 
distributions of atmospheric muons on charm production. An example of 
a model for charm production that is near the upper limit is the RQPM 
model of Bugaev et al. [1998]. Figure [2] shows the angular distribution in 
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Angular dependence of 10 - 20 TeV neutrinos 
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Figure 2: The curved line shows the angular dependence of from decay 
of pions and kaons while the horizontal line shows the isotropic prompt 
component in the model of Bugaev et al. [1998] for E u = 10 to 20 TeV. 

the energy range 10 < E v < 20 TeV using this model. The prompt and 
normal components are shown separately. 

In the model of Bugaev et al. atmospheric neutrinos from charm decay 
become equal in intensity to neutrinos from decay of kaons and pions at ap- 
proximately 100 TeV. A more recent calculation [Enberg, Reno & Sarcevic, 
2008] predicts a level of charm production an order of magnitude lower than 
the model of Bugaev et al. In this case, the crossover occurs about a factor 
of two higher in energy. 

5 Conclusion: detecting neutrinos 

The neutrino effective area is defined so that the product of neutrino inten- 
sity multiplied by the effective area gives the event rate. The probability 
that a muon neutrino on a trajectory that will intercept the detector gives 
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Muon-neutrino effective area for an ideal km 3 detector 
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Figure 3: Muon neutrino effective area for an idealized spherical detector 
with cross-sectional area of 1 km 2 . 

a visible muon in the detector is 

P(E u ,E, Mn ) = N A / d E»R(E„E E ^ in ), (7) 

where the integrand is the product of the charged current differential cross 
section and R is the average distance traveled by a muon with energy En 
at production before its energy is reduced to -E^min, the minimum energy 
required upon entering the detector for the event to be reconstructed well. 
The additional contribution from neutrinos that interact within a large de- 
tector can be computed in a straightforward way and added as an extra 
contribution. 

The neutrino effective area is 

A eS (e,E u ) = e(e)A(e)P(E u ,E flMn )e-^ E ^ NAX(e \ (8) 



s 



where e{9) is the efficiency for a detector of projected area A{9) to detect 
a muon incident at zenith angle 9. The exponential expresses the muon 
attenuation in the Earth for angle 9 below the horizon, where X{9) is the 
amount of matter (g/cm 2 ) along the chord through the Earth. Figure [3] 
shows A e R for an idealized spherical detector with physical projected area 
of one km 2 for several different zenith angles. The shaded region indicates 
the energy range of interest between 10 and 100 TeV. Because the neutrino 
spectrum decreases strongly with increasing energy, most of the rate will 
be near the lower end of this range. In Fig. [2j the prompt contribution is 
about 40% of the normal 7r, K channel contribution to the intensity of 
near cos(#) ~ —0.9, which roughly corresponds to the core region. Because 
the level of charm is uncertain, this source of uncertainty in the angular 
distribution of will be a significant limitation to the measurement of the 
core-mantle transition. Alternatively, a careful measurement of the angu- 
lar distribution of atmospheric neutrino induced muons with good energy 
resolution in this energy range may clarify the level of charm production. 
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